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Abstract: We construct a model for delayed electroweak symmetry breaking that takes 
5^ \ place in a cold Universe with T <C 100 GeV and which proceeds by a fast quench rather 

than by a conventional, slow, phase transition. This is achieved by coupling the Standard 
Model Higgs to an additional scalar field. We show that the quench transition can be made 
fast enough for successful Cold Electroweak Baryogenesis, while leaving known electroweak 
physics unchanged. 
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1. Introduction 

It is well known that electroweak baryogenesis within the Standard Model is subject to a 
double kill: the electroweak phase transition is a cross-over [|^ so that the requirement for 
non-equilibrium is never fulfilled; and at the electroweak phase transition temperatures, 
the CP-violation in the fermion mass matrix is too small by many orders of magnitude 
(see [Q] for a review). However, it has also been known for some time that the observed 
baryon asymmetry could be explained if the electroweak transition is cold rather than 
hot P, H, [5|, 0] with T ^ 100 GeV, and if the transition itself is a fast quench rather 
than adiabatic P, IC]. Moreover, recent analytical work O and numerical simulations 



have indicated that in such a cold environment. Standard Model CP-violation is much 
larger than at electroweak temperatures and is indeed strong enough to account for the 
baryon asymmetry ||8[. In Cold Electroweak Baryogenesis, most of the baryon asymmetry 
is produced at the initial quench when the Higgs field is rapidly falling down the slope of 
the potential. Baryon production essentially stops after the first oscillation, after which 
the coherent Higgs field starts decaying, thereby reheating the Universe. The reheating 
temperature should be low enough that sphaleron diffusion does not subsequently wipe 
out the baryon number. 

In the early Universe the rate of any global change is necessarily related to the Hubble 
rate, which at the electroweak scale is very small. Hence, in the absence of a first order 
phase transition, the realization of a fast quench does not appear to be possible within the 
Standard Model. But if the Higgs field were coupled to some beyond-the-Standard Model 
fields, the situation might change. Indeed, the Higgs could be just one of many scalar fields 
such as in the low energy limit of string theory. 



The inflaton field could also be an example of such "moduli" , and a low temperature 
state and a fast quench could both follow from single field ("inverted hybrid" -type) low-scale 
inflation [|ll], |l^, ^, 13|. However, this is at the expense of allowing a rather general inflaton 



potential (up to the sixth power in the field), as well as some fine tuning of parameters. 
Presumably this tuning can be alleviated by considering more complicated models (see, for 
instance |l^). 

Alternatively, as in the present paper, one can assume that inflation is decoupled from 
the electroweak phase transition and takes place at some high energy scale. We then 
study the generic conditions for the cold, quenched electroweak phase transition in the 
presence of an extra scalar field coupled to the Standard Model Higgs. In particular, we 
focus on the possibility of triggering the quench through the cosmological expansion, via 
thermal corrections to the effective potential. We will demonstrate both that, under certain 
conditions, the phase transition can be delayed until T ~ 0{1) GeV, and that the transition 
can be fast (in a sense to be specified below). 

We will describe briefiy the problem of the fast quench in the next section. In section 
m we couple the Higgs field to another scalar field and show that a fast quench can be 
achieved. In section H we show how to trigger a quench of the a through finite temperature 
corrections to the mass and the expansion of the Universe. We then briefly consider effects 
such as defect formation, bubble nucleation and loop corrections in section |5[ In section y 
we introduce the criteria to be fulfilled in order that the Standard Model physics remains 
unchanged and section contains our conclusions. 

2. The Higgs quench in the Standard Model 

Consider the standard cosmological scenario, with inflation leading to reheating of the 
Standard Model degrees of freedom, and a reheating temperature T^eh 3> 100 GeV. The 
Higgs potential is (neglecting for the moment interaction with fermions and gauge fields), 

Vi<P) = Vo-^<P^ + ^<P''. (2.1) 

At zero temperature, the symmetry is "broken", in the sense that the field acquires a 
vacuum expectation value v = ^u/vA, which in the Standard Model has the value v = 
246 GeV. At finite temperature, the symmetry is "restored" through a thermal mass (at 
leading order in coupling constants), so that 

vHT) = ^' - 'f^^^\ ml^{T)^KT\ (2.2) 

where k < 1 is a constant^, assumed to include interactions with all the Standard Model 
degrees of freedom. When the right-hand side of (2.2) is negative, the expectation value is 



zero. Since the temperature decreases with the scale factor a as T oc 1/a, we have that the 



^This is the T ^ m result, which we use for illustration. A more refined treatment would include the 
effects of a finite mass (see also section W). 



speed of the quench can be quantified as 



' ^'^'.m-/^')T=T=^' (2-3) 



[m. 



where the subscript q denotes the time of the quench, m^^iTg) = fj.'^- At the electroweak 
scale Vq c^ ^ c^ Tq ^ 100 GeV while the Hubble rate is very small so that the quench 
rate is minuscule with Vq = Hq/ ji ~ 10~^^. Numerical simulations suggest that for Cold 
Electroweak Baryogenesis to reproduce the observed baryon asymmetry, we must require 



m 



t;g>0.1, Tg-lGeV. (2.4) 

Clearly, in the Minimal Standard Model this scheme does not work, and the electroweak 
transition is an equilibrium cross-over. Note also that even ignoring the constraint on Vq, 
Tq ~ 1 GeV would require the thermal correction strength k ~ 10^, since the Higgs mass 
m]j = 2^2 ^ (100 - 200GeV)2 is fixed. 

3. Quench in the presence of an additional scalar 

Consider now a system of two scalar fields; the Higgs, (j), and a Standard Model singlet 
field, a. The potential is chosen to be 

Note that the potential is not renormalizable (o"^), as it turns out that a renormalizable 
potential (a^, o"^) does not lead to a fast enough quench (see also [^, |^). Denoting by 
{va,Vfp) the global minimum, we have the constraints 

- g\l - XavI + X&vt = 0, (3.2) 

ml - g'vl + \^vl = 0, (3.3) 

and to avoid a spurious cosmological constant, we impose that in the minimum, the po- 
tential should vanish. 



V{v^,v^) = Q. (3.4) 



This means that 



a'^vl = ml + X^vl, (3.5) 

X^vt = 12^0 + 2mlvl - X^v% (3.6) 

XqvI = 12Vb + Smlvl (3.7) 

We fix v^ = 246 GeV (electroweak physics), Vq = 100^ GeV^ (to end up in the broken 
electroweak phase after thermalisation) and A^ = m?^/{2vV), mn = 160 GeV (Higgs mass 
allowed by experiment). This leaves the free parameters v^ and m^j,. 
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Figure 1: The Vq as a function of Va and m^. 



Imagine now starting off at (/> = 0, a = +e ^ 1. As long as cr < cr„ 



m. 



is 



enforced. Then upon reaching cjq, the effective mass of (f) flips sign, and (j) will go through 
a spinodal transition and electroweak symmetry breaking. We are interested in the speed 
of the quench, and so we calculate 
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2^di 



[rrij 
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■■m^/g 



(3.8) 



where /i = 100 GeV is just a normalisation scale, in order to compare with | IC ] . For effective 
baryogenesis, we require Vg > 0.1. We easily find 



^g — n 



(3.9) 



Simple energy considerations show that at cr„, the velocity of a is 



a ~ 



'2X4 fm^\ 2A6 f fUd^ 



6 V 5 



(3.10) 



Fig. |l] shows Vq in Va- — m^ space. Only the range Vq > 0.1 is included. Let us for the 
purpose of illustration choose a representative point, Va- = 3000 GeV, m^ = 150 GeV (see 
also Fig. |5|). In that case, we have 

v„ = 3000 GeV, m^ = 150 GeV -^ 

Vq = {).l27, A4 = 3.89xl0"^ Ae = 7.25 x 10^^^ GeV~^ 5^ ^ 3.92 x 10"^. 

(3.11) 

In an expanding Universe, before the quench, the dynamics of the rolling field are given 

by 



a + ZHa = A4cr^ - \&a^ ■ 



(3.12) 



Because of Hubble friction, energy conservation as used in ( 3.10| ) is in principle not strictly 
enforced. We can estimate the effect of the friction, say at a = aq by writing, for our 
example point 



3Ha , 



3000 GeV, m^ = 150 GeV -^ - — 3 |^, = 1.8 x W'^^, (3.13) 



and where we have used that H^ ~ Vq/SM^^. We find in general that this ratio is less than 
10~^^. Hence, the force from the quartic term vastly dominates Hubble friction, which we 
can therefore ignore. 

4. Triggering the a quench 

So far, we have simply assumed that the a field happened to find itself at the initial value 
cr = +e <C 1. In fact, the initial condition is the post-inflationary equilibrium Universe 
at finite temperature T. We will now make the following additions to our model, 1) The 
a has a small negative mass term and 2) it is coupled to a (set of) light fields ^. Their 
purpose is to provide a time-dependent mass for a through thermal corrections. Hence, we 
should add to the potential (pj|), 

y(a,0) ^ V{aA,i) = y{^A) - -fo'' + -cy'^e + -j^e. (4.1) 

In that case, the quantum field equations read, in the Hartree approximation, 

d + -iHd -dla = - {-ml + ^(e^) - 3X^{af + IdXeiaY - g'^{(P^)) a, (4.2) 

4> + 3H4>- dl<^ = - (m^ + ?>\^{^f - g\a^)) 0, (4.3) 

I + 3Hi - dli = - {ml + K{a^)) i. (4.4) 

During this stage of the evolution, {a) = 0, and we are in equilibrium. We are interested 
in the case where we wait until the temperature is very low, in which case the expansion 
of the Universe will eventually be dominated by the vacuum energy Vq. This is simply a 
question of waiting long enough. This also explains the choice of Hubble rate at the end 
of the previous section. 

We also note that^ A^ » (7^ » A4 ~ Aev^, and we will set k ~ A,^. Then we have the 
coupled gap equations in equilibrium 

=i^, nt={^tl--i)-\ (4.5) 

(.^)=/!i±iZH, „J=(,</T_i)-, ,4.6) 

Jk ^k ^ ' 

(a=/^ii^. ni=(e"I/--l)-\ (4.7) 



^Certainly for our example point Va — 3000 GeV, m^ = 150 GeV, 



with 

4 = ,f^^^^^^), u:l = ^-ml + K{e), 4 = ^ml + K{a^). (4.8) 

These equations can in principle be solved numerically, although some care has to be given 
to renormalisation. We will not do so here, but concentrate on the important features of 
the system of equations. 

The first equation is precisely the one we considered in the Higgs-only case in section |2|, 
although here we have not assumed a quadratic temperature dependence of the mass. In 
our approximation, (j) decouples from the rest of the system. We will choose m? such that 
in the range of temperatures we are interested in mi ^ k((7^). Then the ^ is just a free 
massive field at finite temperature providing a time-dependent mass to the a field. 

The middle equation (4.3) tells us, that when ■\/k{^'^) = rua, the a field goes through 
a spinodal transition. This is very slow indeed, since the rate of change of the effective 
mass 

i^ = |(_„e+«(a).«.T^ ,4.9) 

is now very small (as for the case of the Standard Model Higgs in section g). 

However, this does not matter, since at some point the quartic term will take over the 
dynamics, roughly when 



I m'i 



'Tave = ^A^)^J^. (4.10) 

We only need to make sure that this is well before the (/> quench, i.e. choose nia such that. 



J'f « ^. (411) 

V ^^4 g 

which in our example amounts to rriu ^ 15 GeV. This choice also means that the location 
of Va- ( |3.2| ) and the estimate of a ( |3.1C ) is unaffected by the introduction of ma- Finally, 



choosing rricr small enough ensures that the temperature of the Universe is well below 
100 GeV, as required by Cold Electroweak Baryogenesis. 

In principle, all IR modes with \k\ < ma will participate in the spinodal "roll-off". For 
our purpose, we will simply let the quantity Uave play the role of o", when considering the 
subsequent dynamics as in section ^. 

Again, we need to make sure that the dynamics will not be Hubble friction dominated, 
and we write 

t;„ = 3000 GeV, m<i = 150 GeV -^ ^-^| =4xl0-l^ m,, = 1 GeV. (4.12) 

Note that although a spinodal transition is triggered the moment M^ < 0, this transition 
is very slow at first, and involves only the very IR modes^. Asymptotically, of course, 
M^ — )• — m^, and the transition will complete, and even before that the quartic term in 
the potential will have taken over. 

^In the approximation (^^) oc T^, one may solve the field evolution exactly by writing M^{t) = 
ml{e~^"* - 1) [[sf or in the linear approximation M^{t) = -ml{2Ht) ^. 
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Figure 2: Contour plot of the energy of the critical bubble in units of the temperature, 
logjQ(i?crit/T), for our example point w^ = 3000 GeV, m^ — 150 GeV. 



5. Defects, bubbles and loops 

We first note that the a potential has Z2 symmetry, and topological defects (kinks) po- 
tentially form in the transition. However, because the transition is so slow at first, the 
density of kinks will be very small. In fact, the dynamics may be so slow that even be- 
fore M^ < 0, a first order phase transition may be triggered, with nucleation of bubbles. 
In particular, a first order phase transition is possible at finite temperature, even though 
the tree-level/zero-temperature potential does not have a "bump". Whether this hap- 
pens depends on the details of the potential, and an accurate determination of the bubble 



nucleation rate requires careful numerical simulations |17]. 

However, a simple estimate of the nucleation rate can be given in terms of the critical 
bubble 

^crit = ^ ~'r::/ , (5.1) 



Nucleation l-' j m4 ^' 

— ~ i e 



V 



crit 



Ays 



where -Ecrit is the energy of the critical bubble, Vjump is the field expectation value that the 



field jumps to inside the bubble, the "latent heat" is AV = V{a 
the wall tension is given by 



0) - V{a 



V 



jump 



) and 



Sfv 



jump; 



V2{V{v) - y(7;finai)) dv. 



Let us first ignore the Higgs and consider the effective potential 



Via) = Vo + 



kT^ 



m„ 



-a 



A4 4 , ^^6 6 



(5.2) 



(5.3) 



with the minimum ■ugnah ai^d where Vq is fixed so that V^(ufinai) = 0. One can find Ujump 
by the over-shoot/under-shoot method, but we will simply calculate the exponent —^, as 
a function of T and I'jumpj for k = 1. We need to compare the nucleation rate within a 
volume V — H^^ to the Hubble rate H. Since a = becomes the global minimum for 
T > 7 GeV, negligible nucleation amounts to 

^'''" » In (^) ^ 130. (5.4) 



Using again our example point, Fig. ^ shows the contours of logiQ{EcTit/T). In fact, the 
lower right triangle corresponds to "tunneling" that does not pass the bump, and so cor- 
responds to nucleation of a subleading bubble which would contract again. 

When including the Higgs, we should consider the option of tunneling directly into 
the global minimum {va,V(j,) = (3000 GeV, 246 GeV). To estimate the impact on the rate, 
we can compare the potential at the global minimum to the potential at the maximum 
tunneling rate region from Fig. ||, around (cr, </>) = (2250 GeV, GeV), 

y(o,o) -y(225o,o) , , 

^ ^ ^ ' -0.91. (5.5) 



y(0,0)-y(3000,246) 



Hence, there is no dramatic gain in "latent heat" from tunneling directly to the global 
minimum, to compensate for the increase in wall tension from tunneling further in field 



space. We conclude that the constraint (5^) is easily met, and we therefore do not expect 
bubble nucleation to play a role in this model. ^ 

Finally, we emphasize that we consider (|3.1[ ), ( ^.1[ ) to be an effective potential, and in 
that sense, quantum corrections are already included. However, as a check of consistency, 
it is relevant to consider corrections to the a mass at the one-loop level, since our scenario 
is based on m^ being negative (and small). Following [13|, the contribution to the effective 
potential is simply 

i 

where i/ is a renormalisation scale, and mf{a) are the eigenvalues of the mass matrix 
(including ^), at a given value of a (and (p and ^). Away from a = 0, we need to include 
( |5.6| ) for all rrii , and we run into the well-known problem of negative eigenvalues (for small 
a in the a direction; and near ag in the (j) direction), making the effective action complex^. 
Given the correct effective potential, we could solve for a different set of "bare" parameters 
to satisfy the original constraints. But at one loop <\5.(\) , the Higgs quench rate constraint 
is beyond our reach. 



^This turns out to be a result of having a rather small A4 and correspondingly a very large iijump, 
compared to for instance a similar calculation in the Standard Model Higgs potential (for Higgs masses 
where there is a first order phase transition). 

^The effective action can be calculated by other methods such as a 2PI resummation or full-fledged 
lattice Monte-Carlo, or even by invoking some Maxwell construction, but this is far beyond the scope of 
the present work. 
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Figure 3: The smaller mass eigenvalue to_ as a function of Va and to^. 



We will therefore restrict ourselves to considering the a mass term around a = (j) = 
^ = 0, and since k ^ (7^ ~ A4, we concentrate on the contribution from the ^ field (which 
is real) . We use mf = mi + kcj"^ to find 



1^ 



=i=o 



umi I I nil 



(5,7) 



We will assume a renormalisation scale v > m^, so that the log is negative. But even 
ignoring the log, as long as m^ < -v/327r/K|mo-|, the effective potential still has negative 
curvature at o" = 0. ^ 



6. Electroweak physics 



In order not to conflict with known electroweak physics, we fixed the (j) vacuum expectation 
value to f (^ = 246 GeV. Because of the (f) — a coupling, there is mixing in the vacuum, with 
a mass matrix that reads 



M^ 



2\^vl -2g^VaV^ 

-Ig^VaV^ -2A4W^ + 4^g^4^ 



(6.1) 



This can be diagonalised in a straightforward way, and Fig. ^ shows m-_ as a function 
of Vfj — rriip, corresponding to the eigenvalues ± of the mass matrix. We see that for 
m^ > 182 GeV, the vacuum becomes unstable, ruling out that part of parameter space. 
For our example point, we have 



3000 GeV, m^ = 150 GeV 



m+ = 164 GeV, m_ = 18 GeV, (6.2) 



®The constraint on m^ introduced above, m^ ^ y/Kia-"^) ~ 1 GeV, was only for practical purposes, and 
is not essential. 



with a mixing angle 



6 = tan 



ruj 



2\^vl 



2g^VaV^ 



0.225, 



e = 13°. 



(6.3) 



One may worry that when the mostly-(/> mode m^ is heavier than the mostly-o" mode 
rri-, reheating wih proceed into mostly-cr particles rather than the Standard Model degrees 
of freedom. However, because of the smallness of g^, this does not happen. For instance, 
the (f) ^ 2a decay mediated by the term 2g'^V(f 



a^ , is 



SirrriH 



4M2 



M 



,t2 



2.26 X lO^^GeV 



(6.4) 



Comparing this to the total width of Higgs decay into Standard Model degrees of freedom 
in this mass range, T ~ (10^'^ — 1) GeV |jl^, we conclude that the vast majority of the 
available energy will be channeled into these. 
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Figure 4: The mostly-c^ to two mostly-a decay rate F as a function of Va and m^. 



7. Conclusion 

Using a simple implementation of the Higgs-cr potential, we have argued that the expansion 
of the Universe can be responsible for triggering a fast electroweak symmetry breaking 
transition. At the same time, this transition can be delayed until the temperature of the 
Universe is far below the electroweak scale. 

The model can be summarised in the potential 



vio, ^) = v„- !f .^ - ^.* + ^,» - p<(' + !^0= + ^^^ + 'i.^e + '^e. 



(7.1) 



and we split the evolution into 3 stages: 

''Expanding around the vacuum, a = v^ + Sa, (j) = v^ + i 
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Figure 5: The region in parameter space consistent witli all three constraints presented in Figs. |i|, 
H and ^. The dot is our example point Vg- — 3000 GeV, m^ = 150 GeV. 



• 1) The post inflationary Universe in equilibrium cools under Hubble expansion until 
the a is driven to a tachyonic transition. 

• 2) When (Jave = (c^) ~ m'^/X^ Hubble expansion can ignored and a rolls down its 
potential, triggering a fast tachyonic quench of the Standard Model Higgs, (j). 

• 3) (/) and a settle in the global minimum of the potential, for which the physics is 
constrained by electroweak phenomenology. 

The cost is the addition of an extra scalar field coupling to some other matter sector 
"^" , and with self-couplings a'^ and a^ . This is partly a result of insisting that electroweak 
physics is unaltered, but it also comes from restricting ourselves to a minimal model. More 
complicated cr-potentials can be envisaged, and for instance replacing (o"^, a^) by {a^, a^) 
or {a^, a^) works as well. The input values used here for the Higgs mass {niH = 160 GeV) 
can be relaxed within the range still allowed by experiment, say mn G [115 : 200] GeV. 
Similarly, the height of the potential Vq can be relaxed to ~ 200^ GeV, while still avoiding 
finite temperature electroweak symmetry restoration after reheating. 

The model is very similar to the low-scale inflation model proposed in ||5|, |l3| , although 
since now a is not required to be responsible for inflation and the CMB fluctuations, a 
simpler potential and weaker tuning is allowed. Still, as in |13] there is (p-a coupling and 



mixing, providing an observational signature of this model. However, in contrast to |]13| , 
one mass eigenvalue (m_) is smaller than ?7i//(~ m^), and the Higgs can in principle decay 
into (T particles. This is however suppressed by a small (j) — a coupling. 



11 



We conclude that although the expansion of the Universe is very slow at the electroweak 
scale, it can be amplified in a rather straightforward way to provide a fast electroweak 
transition. This transition can be delayed indefinitely (essentially until T ~ rricj), and so 
also until the Universe is cold enough that the computation [^ of the effective Standard 
Model CP-violation becomes reliable. And so, although far from being the final word on 
the matter, we consider the scenario presented here an interesting option to realise a cold, 
fast electroweak transition. 
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